Few-Body Systems 0, 1-5 (2010) 



Few- 
Body 
ystems 




© by Springer- Verlag 2010 
Printed in Austria 



Two-boson Correlations in Various 
One-dimensional Traps* 



A. Okopiriska **, P. 



Koscik 



Institute of Physics, University of Humanities and Sciences, Swi§tokrzyska 15, 25-406 Kielce, 
Poland 



Abstract. A one-dimensional system of two trapped bosons which interact 
through a contact potential is studied using the optimized configuration in- 
teraction method. The rapid convergence of the method is demonstrated for 
trapping potentials of convex and non-convex shapes. The energy spectra, 
as well as natural orbitals and their occupation numbers are determined in 
function of the inter-boson interaction strength. Entanglement characteristics 
are discussed in dependence on the shape of the confining potential. 

1 Introduction 

Entanglement as a measure of quantum correlations is investigated in the hope 
of better understanding the structure of strongly-coupled many-body systems. 
Recently there is a growing interest in studying few-particle trapped systems, 
since they became accessible in experiments with ultracold gases in optical lat- 
tices and microtraps. The interatomic interaction can be there considered as a 
contact one. By choosing the transverse confinement much stronger than the 
longitudinal one, the quasi-one-dimensional systems with an effective interaction 
gmS(x2 — xi) of an adjustable strength gn> may be experimentally realized p]. 
In the Tonks-Girardeau (TG) limit of gm —> oo the system is solvable for arbi- 
trary trapping potential [2]. Theoretical consideration of such a system evolution 
from weak to strong interactions is thus of interest. 

We discuss entanglement properties for a system of two bosons interacting 
through a contact potential and subject to a confining potential V[x). The di- 
mensionless Schrodinger equation takes a form 
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where the Hamiltonian reads 

H = ~\& 1 ~1& 2 + V{Xl) + V{X2) + 9id6{X2 " Xl) - (2) 

Since the two-boson function is symmetric and may be chosen real, there exists 
an orthonormal real basis {v{\ such that 

<j)(x 1 ,x 2 ) = ^2 hvi{xi)vi{x 2 ), (3) 
l 

where the coefficients ki are real and ^ kf = 1. Therefore 

/oo 
<j)(x,x')vi(x')dx' = kivi(x), (4) 
-oo 

which means that vi are eigenvectors of the two-particle function. It may be 
shown that vi are also eigenvectors of the density matrix, known as natural 
orbitals. Density matrix decomposition is given by p(x,x ) = ^l v l( x ) v l( x )> 
where the occupancies A; = kf. The number of nonzero coefficients k\ and the 
distribution of their values characterize the degree of entanglement. 



2 Optimized Configuration Interaction Method 

The configuration interaction method (CI) consists in choosing the orthogonal 
basis set in the Rayleigh-Ritz (RR) procedure so as to ensure proper symmetry 
under exchange of particles [3]. For the two-boson system, the CI expansion reads 

4>(x 1 ,x 2 ) = y^aijipij(xi,x 2 ), (5) 

where (xi,x 2 \ij) = ^ij{x\, x 2 ) = bij[pi (xi)ipj(x 2 ) + (pj(xi)cpi(x 2 )] with bij = 1/2 
for i=j and bij = l/y/2 for ij^j. Exact diagonalization of the infinite Hamilto- 
nian matrix H nm ij = {nm\H\ij) determines the whole spectrum of the system. 
Truncated matrices [-H]jVxiV allow determination of successive approximations 
to the larger and larger number of states by increasing the order N. We use the 
one-particle basis of the harmonic oscillator eigenfunctions 



<Pi(x) = f ^? j Biiyffix) exp [-Qx 2 /2] . (6) 

Following the optimized RR scheme [4], we adjust the value of the frequency Q 
so as to make stationary the approximate sum of iV bound-state energies, by 
requiring 

5Tr[H] NxN = 

5Q { ' 

Such a way of proceeding has been shown to improve strongly the convergence 
of the RR method [4j [5] . The iVth order calculation provides approximations to 
many eigenstates, which enables a direct determination of natural orbitals by 
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representing them in the same basis ([6]) as v(x) = YlPnPn ( x )- This turns the 
eigenequation into an algebraic problem 

and a nm are determined from diagonalization of [H]nxN- By diagonalization of 
the matrix [A]j\f X ]\f, the approximate coefficients k n may be determined. Due to 
the fact that Yl A 2 nm = 1, their numerical values satisfy ^ k% = 1. 



J, 

a nn for m = n 

(8) 

2~ 1 / 2 a mn for m ^ n 



3 Results 



In the case of harmonic confinement V(x) = mx 2 /2 and the contact interaction, 
the two-particle wave function may be analytically expressed [6]. This allows 
determination of the occupancies A, = kf by discretizing (J3|). The two largest 
occupancies for the ground state are shown in Fig. [1] in function of g\D- The 
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Figure 1. The occupancies Ao and Ai for a harmonically confined two-boson system in function 
of gm, their TG limits are marked by horizontal lines. 



state is non-entangled (Ao = 1) only if the bosons do not interact. The weakly 
entangled "condensed" state with only one orbital significantly occupied is real- 
ized at very weak interactions, g\D ^ 0.1. With increasing guj, the entanglement 
grows, which shows up in the increase of Ai at the cost of Ao- The occupancies 
monotonically approach their TG limits Ag^ ~ 0.7745 and A^ ~ 0.1765. 

Entanglement properties in the case of multi-well potentials are markedly 
different. Using the optimized RR method, we calculated the natural orbital oc- 
cupancies of ground states in double-well potential V2well(^) = 2?a(l — ax2 ) 2 an d 
triple-well potential V3 we n(x) = ^x 2 — ax 4 + ^-x 6 . The potentials have minima 
of the same depth and the maxima of the same height, controlled by the param- 
eter a. The results for a = 0.025 are plotted in Figs. [2] and EJ where the upper 
left presents the shapes of the potentials, and the lower left, the two largest 
occupancies Ao and Ai in function of g\o- For gi£> = 0, the ground state is non- 
entangled, as Aq = 1. With increasing interactions, Aq decreases and Ai grows, 
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Figure 2. Double-well potential (upper left), the occupancies Ao and Ai in function of gm 
(lower left) and two-boson densities (right) for (a) gm = 0, (b) gm = 2.5 ■ 10 -8 , (c) gm = 
5- 1(T 8 , (d) gio = 10" 6 



monotonically approaching the TG limit of non-entangled "fragmented" state, 
Xq G = Xf G = 0.5. The critical value gf D , above which Ao ~ Ai, is much larger 
for the triple- well potential than for the double- well one. The dependence of the 
two-boson density on g\D for the double-well potential is shown on the right of 
Fig. [21 For noninteracting bosons, the probability of both being in different wells 
is the same as being in the same well. With increasing gio, the probability of 
finding the bosons in the same well quickly decreases and above g^ r D the state is 
almost fragmented. In the triple- well case (right of Fig. [3]) the particles live in the 
middle well, only above gf D the probability of finding a particle in an external 
well becomes considerable. In the TG limit of non-entangled "fragmented" state, 
one particle is localized in the middle and the other in one of external wells. 




Figure 3. Same as Fig[2]but for the triple-well potential. The two-boson densities (right) for 
(a) g 1D = 1, (b) g 1D = 1.97, (c) g 1D = 1.985, (d) g 1D = 2.05 



4 Conclusion 

The optimized CI method proves very effective in determining the spectrum and 
the natural orbitals of the two-particle confined systems. 
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